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ABSTRACT 

We present a complete survey of CO (1 0) emission in the Local Group dwarf irregular IC 10. 
The survey, conducted with the BIMA interferometer, covers the stellar disk and a large fraction 
of the extended H I envelope with the sensitivity and resolution necessary to detect individual giant 
molecular clouds (GMCs) at the distance of IC 10 (950 kpc). We find 16 clouds with a total CO 
luminosity of 1 x 10^ K km s~^ pc^, equivalent to 4 x 10^ M0 of molecular gas using the Galactic 
CO-to-H2 conversion factor. Observations with the ARO 12m find that BIMA may resolve out as 
much as 50% of the CO emission, and we estimate the total CO luminosity as ^ 2.2 x 10^ K km s~^ 
pc^ . We measure the properties of 14 GMCs from high resolution OVRO data. These clouds are very 
similar to Galactic GMCs in their sizes, line widths, luminosities, and CO-to-H2 conversion factors 
despite the low metallicity of IC 10 {Z « 1/5Zq). Comparing the BIMA survey to the atomic gas and 
stellar content of IC 10 we find that most of the CO emission is coincident with high surface density 
H I. IC 10 displays a much higher star formation rate per unit molecular (II2) or total (H I + II2) 
gas than most galaxies. This could be a real difference or may be an evolutionary effect — the star 
formation rate may have been higher in the recent past. 

Subject headings: ISM: molecules, galaxies: dwarf, galaxies: ISM, stars: formation 



1. INTRODUCTION 

Star forming dwarf galaxies tend to have low metal- 
licities, intense interstellar radiation fields, and shallow 
potential wells. For these reasons, dwarfs are often used 
as astrophysical laboratories in which to study the ef- 
fect of extreme conditions on the interstellar medium 
(ISM) and star formation. The similarity between lo- 
cal dwarf galaxies and the first star forming systems — 
which were also low mass and chemically primitive — fur- 
ther motivates such studies. Unfortunately, the smallest 
actively star forming systems are not very luminous and 
are therefore difficult to observe at any significant dis- 
tance. Systems the size of the Large Magellanic Cloud 
(LMC) may be studied in some detail out to 10 Mpc or 
more, but detailed observations of smaller systems are 
possible only in the Local Group. In practice this limits 
such studies to a handful of systems, the Local Group ir- 
regular galaxy IC 10, the SmaU Magellanic Cloud (SMC), 
NGC 6822, and perhaps a few other Local Group dwarfs 
(jkatco I99J). 

In this paper we present a new study of the molecular 
gas component of IC 10. We have conducted a com- 
plete survey of ^^CO (1 0) emission from IC 10 using 
the BIMA interferometer. This survey covers the optical 
disk and much of the extended H I structure surround- 
ing IC 10 with the resolution and sensitivity necessary 
to detect individual giant molecular clouds (GMCs). We 
also present new single dish observations towards most of 
the CO emission from IC 10 using the ARO 12m and we 
measure macroscopic properties of 14 GMCs us i ng th e 
high-resolution OVRO observations by IWalter I 1)200,*^ . 
We combine these data with literature observations of 
IC 10 at several wavelengths to address three questions: 
1) How do the GMCs in IC 10 compare to the GMCs 
found in Local Group spirals? 2) Can we predict the 
molecular gas content from the hydrostatic pressure? 3) 



Do stars form out of molecular gas at the same rate in 
IC 10 as in spiral galaxies? 

There are several reasons to expect the relationship 
between atomic gas, molecular gas, and star formation 
in IC 10 might be different from that in spiral galax- 
ies. Dust plays a crucial role in setting the abundance 
of molecular gas by shielding molecular gas from dis- 
sociating radiation and serving as the site of molecular 
hydrogen formation. With its low metallicity {Zjc lo ~ 
1/5 Zq, Garnctt 1990), IC 10 might be expected to have 
to have a low dust-to-gas ratio. The interstellar radia- 
tion field (ISRF) in IC 10 should also be more intense 
than in most spiral galaxies due to the high star forma- 
tion rate, low dust abundance (lower extinction), and low 
metallicity (less line blanketing). A more intense ISRF 
will dissociate molecular hydrogen faster and perhaps in- 
hibit the formation of H2 by heating the grain surfaces. 
It may also affect different molecular species differently 
as, for example, H2 self-shields while CO is dissociated. 
Other factors, for instance the lack of shear (perhaps 
diminishing rotational support against CMC collapse), 
the absence of spiral density waves (the sites of CMC 
formation in spirals), and the simple underabundance of 
carbon and oxygen (and thus, presumably, of CO to trace 
and cool molecular gas), further distinguish IC 10 from 
the Galactic environment. 

This paper is organized in the following way. In §2 
we summarize the macroscopic properties of IC 10 and 
present maps of several components of the galaxy. We 
emphasize the similarity of IC 10 to the SMC in gross 
properties and discuss IC lO's vigorous present day star 
formation. In §3 we present our BIMA and ARO 12m ob- 
servations, discuss how we identified signal in the survey, 
describe our algorithm for measuring CMC properties, 
and note several other datasets used in this paper. In §4 
we present the results of the BIMA survey and our CMC 
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property measurements. In §5 we present our analysis 
of the molecular content in IC 10. We emphasize the 
surprising similarity between IC 10 clouds and GMCs 
in large spiral galaxies. We examine quantitative rela- 
tionships between atomic gas, molecular gas, and star 
formation. We consider the hypothesis that the hydro- 
static pressure predicts the molecular to atomic gas ratio 
and look at the efficiency with which molecular gas forms 
stars in IC 10. In §6, we summarize our findings and sug- 
gest some future avenues of investigation in IC 10. 

2. DESCRIPTION OF IC 10 

Along with the SMC, IC 10 is the best example of 
a low-mass, metal-poor, actively star-forming galaxy in 
the Local Group. Table ^summarizes a number of prop- 
erties of IC 10. IC 10 has stellar, atomic gas, and dy- 
namical masses of ~ 4, 2, and 15 x 10^ Mf^, respec- 
tively fJa rrett et all 120031 iHlIchtmeier fc RichteJ ITosl 
[Mateo 199^; by comparison, the SMC has stellar, atomic 
gas, and dynamical masses of ^ 4, 4 , and 24 x 10^ M© 
(jStanimirovic. Stavelev-Smith. fc Jones 2004.). IC 10 has 
a me t allicity of 12 - f log 0/H ft! 8.25 (jLeaueux et al.l 
119791 iGarnetil 1199(1 . intermediate bet ween the SMC 
and the LMC (8.0 and 8.4, respectively. lDufoud[T98l . 
Like the SMC, IC 10 has irregular morphology, ongoing 
high mass star formation, and an extended H I envelope. 
IC 10 is probably associated with M 31 but the sepa- 
ration between the two galaxies is fairly large. At our 
adopted distance (950 kpc), IC 10 is separated from M 31 
galaxy by > 250 kpc. Another close cousin to IC 10 exists 
just beyond the Local Group — the post starburst dwarf 
NGC 1569 at a distance of ~ 2 Mpc. That galaxy is 
also relatively isolated, with atomic gas and total masses 
of ~ 1 and ~ 3 X lO'* (Isr aeilll988D. respecti vely, and a 
metallicity of ~ 1/4 Zq iCalzetti et al.l ll994l 

IC 10 has a higher star formation rate (SFR) than the 
SMC, though the exact SFR is somewhat uncertain. The 
Ha flux implies a SFR of ~ 0.2 Mq yr"-"^ with a stan- 
dard correction for internal extinction, but this value 
may be as high 0.6 Mq yr~^ if the larger extinction 
estimates in the literature are correct l|^ il dc Paz et al. 
l2003tlYang fc Skillmanll l993: B orissova et alj 2000'1. By 
comparison. iWilke et al.i l|200JO estimates an SFR of 
~ 0.05 Mq yr-i in the SMC (with the FIR and Ha 
in agreement after corrections for extinction and absorp- 
tion). Masscy fc Holmes ( 2002), Crowthcr ct al. ( 2003), 
and others have noted IC lO's prodigious content of 
Wolf Rayet (WR) stars. More than 30 have been 
spectroscopically confirm ed ijCrowther et alJ l2003() and 
iMassev fc Holme^ l|2002|) estimate from photometry that 
the total number may be ~ 100 . By contrast, th e 
SMC contains only 12 WR stars l|Massey et aLll2003t) . 
The number of WR stars in a galaxy should provide a 
distance independent estimate of the high mass stellar 
content (and thus formation rate). Thus Ha and WR 
star counts suggest IC 10 to have a star formation rate 
3-4 times th at of the SM C. Measurements of inter- 
nal e xtinction l|Yang fc Skillman .1993.: .Borissova et al.l 
200011 and the higher nu mber of WR stars estimated by 
Massev fc Holmea l|2002|) suggest the SFR may easily be 
as high as ^ 0.5 Mq yi~^. We note that the FIR and 
radio continuum yield SFR e stimates conside rably lower 
(~ 05 Mp yr- ^ lThronson et al...l990; White fc Becked 
Il99l lBeIill2003l) perhaps as a result of a higher UV es- 



cape fraction or lower dust abundance. 

The unusually high content of WR stars has led to 
claims that IC 10 is the nearest starburst galaxy — the 
surface density of WR stars is the highest of any local 
group galaxy and the average across the entire galaxy 
is compar able to the mo st active ly sta r forming regions 
in M 33 l|Massev fc Arm androfi 11995(1 . The disrupted 
morphology of the H I distribution also suggests a vio- 
lent recent history. The H I distribution across the disk 
of the galaxy is characterized by seven large holes, possi- 
bly c arved out by winds or supernovae ifWilc ots fc MilleJ 
[1993). The central regions of these holes are free of 
H I emission down the to a sensitivity limit of 1 x 10^^ 
cm^^ (3cr), equivalent to about 0.1 Mq pc^^. IC 10 has 
a substantial molecular gas content compared to other 
galaxies of its size, suggesting that the period of vigor- 
ous star formation may be ongoing. IC lO's CO lumi- 
nosity of ~ 2 X 10^ K km s~^ pc^ is an order of mag- 
nitude higher than that of the SMC (~ 8 x 10'' K km 
s-i pc2,^izuno ct al."2001) or NGC 1569 (- 10^ K km 
s^^j_£C^jGreyc ct al. 1996; Taylor ct al. 1999) ■ Indeed, 
iWilcots fcM illcr (1998) calculated only ~ 7 - 21 super- 
novae would be needed to create all of dramatic holes in 
the H I distribution; thus, based on the observed num- 
ber of WR stars — possibly as high as 100 and each 
representing a future supernova — IC 10 may have only 
experienced a small fraction of the Type II supernovae 
in store for it in the near future. This would seem to 
place IC 10 in contrast to the comparatively quiescent 
SMC and the post starburst NGC 1569 (where a burst 
of s tar form ation may have ended as recently at 5 Myr 
ago lGreggio et al. 1998). 

Figures ^ and [2| show maps of the atomic gas con- 
tent, stellar surface density, and Ha surface brightness 
in IC 10, with the positions of spectroscopically con- 
firmed WR stars (Crowthcr ct al. 2003) noted in Fig- 
ure m (the larger , ^ 1 00, number of WR stars from 
IMassev fc Holme^ 120021 await spectroscopic confirma- 
tion and publication of thei r locations). The comp lex 
H I structure is described bv IWilcots fc Milleil ^0) in 
the following way: the stellar disk lies within a larger H I 
disk which shows rotation aligned with the stellar disk, 
the whole galaxy lies within a much more extended H I 
structure that is counter-rotating and complex. The H I 
disk extends beyond the stars to the east of the galaxy 
and forms part of a contiguous position-velocity struc- 
ture with the extended H I envelope to the south and 
east of the galaxy. This large envelope may be have re- 
cently interacted with the galaxy in a way that triggered 
the present star formation in the galaxy. Most of the 
star formation lies in the central part of the disk, where 
a large H I cloud is the site of much of the ongoing star 
formation activity. The dramatic holes that give the disk 
its disturbed (almost spiral) morphology are probably a 
result of stellar winds or perhaps supernovae. 

Because it lies close to the Galactic plane (6 = —3.3), 
the distance to IC 10 remains uncertain. Current esti- 
mate s of the distance range from a lower limit o f 500 
kpc l|Sakai et al.l[l999|) to 950 kpc (|Hunteil[200l[) . In 
this paper, we adopt 950 kpc, the distance obtained by 
Hunter ( 2001) using the tip of the red giant branch and 
a reddening confirmed by comparison with color magni- 
tude diagrams. The resulting distance is close to 1 Mpc, 
a value frequently adopted in the literature, making for 
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easy comparison. The values presented in Table ^ and 
the results of this paper have been scaled to this distance 
of 950 kpc. 

3. OBSERVATIONS 

In this paper we present two new datasets: a complete 
survey of the disk of IC 1 obtained with th e BIMA in- 
terferometer (described bv lWelch et alJl996j) . and single- 
dish pointed observations obtained with the Arizona Ra- 
dio Observatory (ARO) 12m to check for emission missed 
by the BIMA survey. In this section we describe the ac- 
quisition and reduction of both data sets and the algo- 
rithm we used to identify signal in the BIMA survey. We 
also introduce several previously published data sets that 
we use in our analysis and summarize our algorithm for 
measuring GMC properties. 

3.1. The BIMA Survey 

The BIMA survey was conducted in the most compact 
BIMA configuration, the D array, and has a resolution 
of 14". The survey took place over the course of four 
observing seasons: fall 2000, spring 2001, fall 2001, and 
spring 2002. We observed IC 10 during 34 tracks rang- 
ing from 2 to 13 hours in length for a total of 250 hours 
of observing time. For most tracks, we also observed a 
planet once or twice to check the absolute flux of the 
phase calibrator. Each track consisted of ~ 20 fields ar- 
ranged on a hexagonal grid with pointing centers spaced 
by 78" (BIMA has a « 100" half-power field of view at 
115.27 GHz). The correlator configuration varied some- 
what over the survey (the central velocity and channel 
width changed slightly) but a typical setup covered a 
bandwidth of 100 MHz (w 260 km s'^) near the velocity 
of IC 10 with a velocity resolution of « 1 km s^^. The 
final survey has a resolution of 3 km across a 150 
km s~^ bandwidth near the velocity of IC 10, so these 
variations in the correlator setup do not affect the final 
data. For the survey 1 K = 2.04 Jy beam 

We reduced the observations using the MIRIAD soft- 
ware package^ . We corrected the observations of the 
phase calibrator and the source for line length variations 
(BIMA monitors the electronic path length to each an- 
tenna by periodically sending a signal to the antenna and 
back and measuring the phase on return, allowing vari- 
ations to be removed from the data during reduction). 
We flagged data with shadowing or very high amplitudes, 
and channels at the edge of the correlator. We adopted 
the flux for the phase calibrator shown in Table 12 and 
self-calibrated on it, assuming that the phase calibrator 
was unresolved by our 14" beam (we see no evidence of 
extended structure in our data). We transfered the gains 
and phases as a function of time to the source. 

We combined the (u, v) data from all 34 tracks, applied 
a natural weighting scheme, and inverted them into a 
spectral line map with velocity channels 3 km wide. 
We applied a CLEAN algorithm to each plane of the 
data cube in order to remove artifacts generated by in- 
complete {u,v) coverage. We capped the algorithm at 
500 iterations in each map and did not clean sources less 
than 2(7 in peak brightness. The final maps cover an 75 
square arcminutes at better than 0.2 K sensitivity with 
an angular resolution of 14", a velocity resolution of 3 
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km s ^ , and a velocity coverage spanning LSR velocities 
from -400 km s^^ to -250 km s~^ 

The flux of our phase calibrator, 0102-1-584, varied by 
a factor of 3 over the two year duration of the survey, 
and as a result, the amplitude calibration of the survey 
may be somewhat uncertain. We used values interpo- 
lated from the BIMA calibrator monitoring campaign, 
adjusted by as much as 30% based on comparisons be- 
tween 0102-f584 and planets in our own observations. 
Based on the variation of the phase calibrator and the 
scatter within the planet/phase calibrator fluxes within 
our own data, we estimate that the flux calibration of 
the survey is uncertain by 30%. These gain errors may 
be nonuniform across the survey because different tracks 
contribute to different parts of the map. 

3.1.1. Sensitivity of the Survey 

On average, 50% of the time in a given track was spent 
integrating on IC 10. The remainder of the time was 
spent on calibration and slewing. The 250 hours of ob- 
serving time used for the survey translates into « 125 
hours of integration on source. The total area targeted 
was 13' X 8' « 100 square arcminutes. This corresponds 
to an integration time per pointing of 2 hours. The 
theoretical Icr noise for 2 hours of observation with the 
D array at 3 km s~^ velocity resolution and our typi- 
cal system temperature (500 K) is ^ 0.05 K. Figure 01 
shows a map of the sensitivity of the survey, measured 
from the RMS variation in the signal-free channels at 
the edges of the map, which is somewhat higher than 
the estimate above due to missing antennas, atmospheric 
decorrelation, and other inefficiencies. About 40 square 
arcminutes of our survey have a la sensitivity better than 
0.1 K, another 35 square arcminutes have Icr sensitivities 
between 0.1 and 0.2 K. Thus the final survey covers 75 
square arcminutes at better than 0.2 K sensitivity. The 
noise in the survey is quite Gaussian. We find 15.9%, 
2.3%, and 0.13% of the emission in the final cube to lie 
below —1(7, — 2cr, and — 3cr, respectively, almost exactly 
what would be expected from a normal distribution. 

3.1.2. Signal Identification in the Survey 

We identified signal using the statistic of merit equal 
to the the product of the probability of generating the 
observed signal or higher in each of five adjacent velocity 
channels along a line of sight (referred to by the sub- 
scripts 1 to 5) from Gaussian noise. This statistic, X, is 
calculated via the formula 

X ^ P{ — ) X P{ — ) X P{^) X P{ — ) X P{^) (1) 

(7 (7 (7 (7 (7 

where P{^) is the probability of generating an intensity 
equal to or greater than li from a distribution of nor- 
mally distributed noise with standard deviation cr. For 
Gaussian noise, P{^) is given by 

I 1 {I^< 0) 

The factor of ^ results from normalizing the error func- 
tion so that erf (oo) — 1 (rather than i). We consider all 
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TABLE 1 
Properties of IC 10 



Property 



Value 



Reference 



Hubble Type 

Distance 

Axis Ratio 

Dynamical Mass 

Absolute B Magnitude 

Stellar Mass 

H I Mass 

Ha Luminosity 

FIR Luminosity 

1.4 GHz Luminosity 

Metallicity 



dirr 
950 kpc 
0.67 

~ 1.7 X 10^ M0 
-16.3 

4 X 10* Mq 
2 X 10* Mq 

1 X 10*° erg s-^ ''•^ 

5 X lO'' Lq 

3.4 X IQi^ W Hz-i ^ 
log O/H + 12 = 8.25 



Mateo QL998) 
Hunter (2001) 
Tarr ett ot ah (2003) 
Mateo C 1998) 
Gil de Paz et al. (2003') 
Jarrett et al. (2003) 
Huchtmeier fc Richtcr (1988) 
^il de Paz ct al. (20031: Thronson et al. ^990) 
Mclissc & Israel ( 1994) 
White & Becker ('1992) 

Teaucux ct al. (1979): Garnett (1990) 



SFRhc 
SFRi.4 Qfj^ 



0.2 Mq yr^"""^ 
0.05 Mq yr-l 
0.07 Mq yr-l 
0.1 Mq yr~^ kpc" 



Gil de Paz ct al. (2003); Kcnnicutt et aL (1994) 
Mclissc fc Israel (1994) ; Boll (20Q3.) 
White fc Becker (1992): Bell (2003) 



this work ( i|4.2i 



CO Luminosity 
Best Estimate 
Lower Limit 
Upper Limit 



2.2 X 10^ K km s"! pc^ 
1.6 X 10^ K km s-i pc^ 
2.8 X 10^ K km s"! pc^ 



''Value has been adjusted to our assumed distance of 950 kpc. 
''Quantity adjusted to use reddening of E(B — V) = 0.77. 
'^Includes 1.1 magnitudes of internal extinction. 
''Ha SFR divided by optical size. 



TABLE 2 
The BIMA Survey 



Observing Season 


Flux of Calibrator 


Dates of Tracks (Duration in Hours) 


Fall 2000 (45h) 


2.5 Jy 


11 Oct (10)=', 12 Oct (12)=^, 16 Oct (10), 22 Oct (13) 


Spring 2001 (75h) 


1.3 Jy 


7 Jun (10), 10 Jun (10), 12 Jun (11), 16 Jun (10), 
17 Jun (10), 21 Jun (5), 23 Jun (4), 29 Jun (4)'', 
1 Jul (3), 3 Jul (8)=^ 


Fall 2001 (83h) 


2.0 Jy 


23 Sep (7), 26 Sep (13), 28 Sep (6), 2 Oct (6), 
3 Oct (11), 6 Oct (5), 7 Oct (5), 12 Oct (8), 
13 Oct (5), 20 Oct (8), 21 Oct (9) 


Spring 2002 (71h) 


3.0 Jy 


25 May (8), 26 May (10), 27 May (5), 28 May (8)=", 
29 May (3), 31 May (8), 6 Jun (9), 8 Jun (2), 
13 Jun (7), 15 Jun (5), 16 Jun (6) 



''Data not used. 

bTrack used different calibrator (0228+673). 



negative intensities to be the result of noise and so assign 
those data a P{^) of 1. Therefore negative intensities 
can never contribute to a detection (which are identified 
by their low values of X). 

We used this statistic to identify lines of sight with 
significant emission in the BIMA survey. We calculated 
X for all combinations of 5 adjacent velocity channels in 
the survey. We then select all regions with X < 10^^ 
and RMS sensitivity of 0.2 K or better. This value of 
X corresponds to « 2a emission across 5 channels, a 
single channel containing « 6cr emission, or a range of 
intermediate cases. We chose this value, X < 10~^, so 
that we do not expect a false detection over the BIMA 
survey. We checked our expected false positive rate by 
applying this algorithm to the negative part of the data 



set (which should consist only of noise). The algorithm 
identifies no significant emission in the negative part of 
the data cube. 

We constructed a mask from the signal we identified. 
The mask is a set of flags (ones and zeros) that identifies 
regions of significant emission in the data cube (a 1 for a 
region containing emission, a for a signal free region). 
In the analysis below we use the emission measured in 
the BIMA survey multiplied by the mask, so that signal 
free regions are treated as having zero emission while 
the data is used in regions containing significant signal. 
To be included in the mask, we required that region of 
the survey to have X < 10~^ across most of a beam 
(at least three adjacent quarter-beam-sized pixels). In 
order to ensure that we did not miss extended, low-lying 
emission around the CO peaks, we included all regions 
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Fig . 1. — H I emission from IC 10 (grayscale, 'Wilcots & Miller 
119981) with stellar surface density contours overlaid (black, derived 
from X-band data from ijarrett et al. 2003|). The stellar surface 
density contours are at 50, 100, 200, and 300 Mq pc-^. The H I 
distribution extends well beyond the optical disk, stretching for 
more than 20 arcminutes (> 5.5 kpc). H I contours show surface 
densities of 1, 10, 20, and 30 Mq pc~^. The surface densities shown 
are line of sight surface densities, with no correction for inclination. 
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Fig. 2.— H« emission from IC 10 IGil de Paz et al."2003') with 
H I contours (10, 20, and 30 Mq pc-^, ^ilcots & Miller 1998) 
overplotted. The locat ions of spectroscopically confirmed WR stars 
ICrowther et al.ll20(13l are shown for comparison. 



within three velocity channels (a typical cloud line width) 
and one and a half beams of each significant region. As 
a result, the mask contains some noise at the edge but it 
is less likely to miss extended emission. The integrated 
intensity from the survey multiplied by the mask is shown 
in Figure^ 

This method does not produce a single point source 




Fig. 3. — Sensitivity of the BIMA D-array survey in Kelvin. The 
wh ite contours show H I surface density (10, 20, and 30 Mq pc~^, 
IWi lcots & Miller 1998). The black contours show sensitivities of 
0.3 and 0.6 K km s^^ per 3 km s~^ channel. 



sensitivity, but in worst case — 5 channels at 2a sig- 
nificance in the region of 0.2 K sensitivity — we find 
all emission with an integrated intensity above 6 K km 
s^^. This corresponds to 1.4 x 10^ over a BIMA syn- 
thesized beam at a Galactic CO-to-H2 conversion factor 
(appropriate for the IC 10 GMCs but perhaps not for 
the diffuse gas). We detect emission with narrow line 
widths in high sensitivity regions down to even smaller 
luminosities, corresponding to masses of < 10^ Mq (1.8 
K km s~^, or 4 x 10^ Mq over a BIMA beam, in the 
best case). A typical sensitivity over the whole survey is 
therefore 10^ Mq. 

3.2. ARO 12m Observations 

We observed 22 pointings in IC 10 with the Arizona 
Radio Observatory (ARO) 12m telescope at Kitt Peak. 
This telescope has a 55" (~ 250 pc) half-power beam- 
width at 115.27 GHz. The data were acquired during 
an observing run in 2002 May 8-16. The locations of 
the ARO 12m pointings are listed in Table |3| and plotted 
along with the results of the BIMA survey and the H I 
contours in Figure — though note that four of the 
pointings lie outside the map. The ARO 12m pointings 
were selected for two reasons. First, they were intended 
to check the flux measurements and line widths of the 
BIMA survey — especially to search for the presence 
of extended emission that might be resolved out by the 
interferometer. Second, they were selected to provide an 
independent search for CO emission at several points of 
interest not detected by BIMA — including H I peaks in 
the extended emission (possible locations of CO beyond 
the optical disk) and tw o of the H I holes identified by 

IWilcots fc MHIea (pgi). 

We observed both polarizations either with the 1 MHz 
filter banks or with the millimeter autocorrclator. We 
observed most pointings for about one hour, though sev- 
eral pointings were bright enough to be detected in only 
~ 30 minutes. We position switched every six minutes, 
with the reference position separated from the pointing 
by 3 arcminutes in azimuth. In only one case, pointing 
14, do we find evidence of emission in the reference posi- 
tions close to the velocity of IC 10. Every six hours, after 
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Fig. 4. — The integrated intensity map resulting from the BIMA D array survey (grayscale). Black contours show the 1, 10, 20, and 30 K 
km surface brightnes s contours. Numbers in dicate the corresponding entry in Table El Light contours indicate H I col umn densities of 
10, 20, and 30 Mq pc~^ IWilcots fc Millerll998fl and the locations of spectroscopically confirmed WR stars are indicated ICrowther et al] 
CO emission is almost exclusively found along lines of sight with Sji/j > 10 Mq pc~^ (.^(H I) = 1.25 X 10^^ cm~'^). 



sunset, and after sunrise, a planet or other strong contin- 
uum source was observed to optimize the pointing and 
focus of the telescope. The median system temperature 
was 345 K. 

We reduced the spectrum for each six-minute scan in 
the following manner. We removed noise spikes and bad 
channels by flagging all channels in each six minute scan 
with absolute values above the 5cr level (none of our 
sources were nearly this bright in a single scan). Several 
channels were known to be bad a priori and we flagged 
these as well. We then subtracted a linear baseline from 
the spectrum and binned it to a resolution of ~ 4 km s~^. 
Finally, we averaged both polarizations and all scans to 
produce the final spectrum for each position. 

The results of the ARO 12m pointings are summarized 
in Table El For spectra with detected signal we quote 
the peak temperature and the central velocity, velocity 
width, and integrated intensity from a three parameter 
Gaussian fit. For spectra without detections we quote 
a 3(7 upper limit to the integrated intensity assuming a 
source 15 km s~^ wide (a typical FWHM line width for 
our detections). For spectra with detected emission, we 
also quote the results from the BIMA and OVRO data 
convolved to the resolution of the ARO 12m. We discuss 



the results of the comparison among the three datasets in 
Section im The final column indicates the BIMA survey 
clouds associated with each pointing. 

3.3. Other Data Used in This Paper 

3.3.1. OVRO High Resolution CO Maps 

We use the high reso l ution OVRO observations of IC 10 
presented bv iWalteTI l)2003j) to study spatially resolved 
GMCs and as a third point of comparison for CO fluxes 
and line widths. These data have an angular resolution 
(FWHM) of 3.4" (13.5 pc) or 4.9" (19.5 pc), depending 
on the location in the galaxy, and a velocity resolution of 
« 0.65 km s~^. The high resolution and good signal-to- 
noise ratio in these data allow us to measure the prop- 
erties of individual GMCs. The OVRO primary beam is 
~ 50" at 115 GHz, so these data are not sensitive at all 
to spatial scales more extended than ^ 25". 

3.3.2. Ha Imaging 

We use an Ha image of IC 10 from iGil de Paz et alJ 
|2003) to measure the star formation per unit area 
(SsFJ?) across IC 10. We correct the Ha intensity for 
the effects of exti nction and then c o nvert to SsFi? using 
the calibration bv lKennicutt et aP l)1994fl . 
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TABLE 3 
ARO 12M Results 



Pointing 


Telescope 


a 


5 


1(7 Noise^ 




Vctr 






Clouds 






(J2000) 


(J2000) 


(mK) 


(K) 


(km s~i) 


(km s"-"-) 


(K km s~l) 


Overlapped 


1 


ARO 12m 


00 20 46.4 


59 18 59.9 


20 


0.10 


-335.7 


8.2 


1.7 


B16 




BIMA"^ 








0.04 


—332.4 


4.2 


0.4 






OVRC^ 








0.03 


—332.0 


3.0 


0.2 




2 


ARO 12m 


00 20 22.3 


59 21 16.9 


19 


0.26 


-329.6 


3.9 


2.6 


B9 


2 


BIMA"^ 








0.13 


—330.0 


3.9 


1.3 




2 


OVRO"^ 








0.05 


—330.0 


4.2 


5 




3 


ARO 12m 


00 20 19.4 


59 20 50.4 


26 


0.20 


-332.6 


5.0 


2.5 B6 




3 


BIMA'^ 








0.08 


—337.4 


4.6 


0.9 




3 


OVRO*^ 








0.04 


—336.7 


3.6 


0.4 




4 


ARO 12m 


00 20 16.4 


59 18 45.5 


14 


0.06 


-332.6 


7.0 


1.0 


B5 


4 


BIMA'^ 








0.04 


—337.1 


3.8 


0.4 




4 


OVRO'^ 








0.04 


—333.9 


2.6 


0.3 




5 


ARO 12m 


00 20 21.3 


59 17 10.6 


36 


0.32 


-337.8 


4.5 


3.3 


B8 


O 


BIMA*^ 








14 


— ooo. / 




1 2 




i) 


OVRO"^ 








0.23 


—339.4 


3.3 


2.2 




6 


ARO 12m 


00 20 27.9 


59 17 01.0 


31 


0.31 


-330.5 


11.1 


8.1 


Bll, B10'= 


O 


BIMA^ 








u.zu 


— oou.o 




9 Q 




g 


OVRO"^ 








0.20 


—330.1 


3.3 


3.3 




7 


ARO 12m 


00 19 54.2 


59 17 22.6 


9 


0.08 


-372.8 


2.4 


0.6 






BIMA'^ 








0.06 


—373.9 


3.7 


0.6 




7 


OVRO'^ 








0.07 


—371.9 


1.6 


0.3 




8 


ARO 12m 


00 19 58.6 


59 18 43.2 


24 


0.10 


-366.4 


3.6 


1.1 


Bl, B2° 


g 


BIMA"^ 








0.06 


—368.4 


3.7 


0.6 




g 


OVRO"^ 








0.03 


—366.5 


2.6 


0.2 




9 


ARO 12m 


00 20 00.3 


59 19 05.9 


15 


0.12 


-368.5 


4.0 


1.1 


B2, Br 


9 


BIMA'^ 








0.07 


-368.5 


3.9 


0.7 




9 


OVRO'^ 








0.03 


—366.2 


2.0 


0.2 




10 


ARO 12m 


00 20 11.3 


59 20 01.3 


17 


0.08 


-343.1 


3.4 


0.7 


B4 


10 


BIMA'^ 








0.07 


-346.0 


3.6 


0.6 




10 


OVRO'= 








0.03 


-343.3 


1.9 


0.2 




11 


ARO 12m 


00 20 39.2 


59 21 12.1 


19 








< 0.51 




12 


ARO 12m 


00 20 26.5 


59 19 16.7 


18 








< 0.49 




13 


ARO 12m 


00 20 19.4 


59 17 57.5 


17 








< 0.46 




14 


ARO 12m 


00 20 00.8 


59 16 41.8 


10 








< 0.27^^ 


B3 


14 


BIMA<= 








0.07 


-364.7 


1.8 


0.5 




14 


OVRO<= 


















15 


ARO 12m 


00 19 55.6 


59 18 01.1 


11 








< 0.30 




16 


ARO 12m 


00 20 11.4 


59 17 01.5 


12 








< 0.32 




17 


ARO 12m 


00 20 40.9 


59 17 41.9 


14 








< 0.38 




18 


ARO 12m 


00 20 30.0 


59 13 00.0 


20 








< 0.54 




19 


ARO 12m 


00 20 35.0 


59 09 42.4 


21 








< 0.57 




20 


ARO 12m 


00 20 19.4 


59 08 35.5 


13 








< 0.35 




21 


ARO 12m 


00 19 21.2 


59 21 13.8 


15 








< 0.40 




22 


ARO 12m 


00 20 51.5 


59 20 19.3 


13 








< 0.35 





^In a 5-2 km s"^ channel. For spectra with a different channel width the noise is adjusted to this channel width. 

^Upper limits arc 3(7 upper limits for a source with velocity width 15 km s ^. 

^All BIMA and OVRO data have been convolved to match the resolution of the ARO 12m. 

■^Emission in the off position. 

*^Cloud not wholly within the beam but may contribute to the ARO pointing. 



SFR 



Mr. 



L{Ha) 



(3) 



yr J 1.26 x lO^i erg s"! ' 

which assumes a Salpeter IMF (from 0.1 to 100 Mq). 

Because IC 10 lies near the Galactic plane, foreground 
extinction is important. We assume a reddening to- 



wards IC 10 of E{B -V) = 0.77 (iM assev fc Armandrofil 
[1995: Hunter 20Q1) . which corresponds to an i?-band 
extinctio n of 2 magnitud es for a Galactic extinc- 
tion law ijCardelli et al.llT989j) . This extinction is con- 
sistent with the Galactic dust c olumn near IC 10 in 
each direction in the dust map of iSchlegel et alJ l[T? 
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(K km s-') 
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Fig. 5. — The integrated intensity map resulting from the 
BIMA D-array survey and the locations of the UASO 12m point- 
ings. The ARO 12m pointings cover most of the CO emission from 
IC 10. Light contours indicate J^hi = 10, 20, and 30 Mq pc~^ 
JWilcots & Miller 1998). Several ARO 12m pointings are outside 
this figure, and appear along the edges. The extent of the survey 
is indicated by the solid black lines. 



Our value is lower t han the extinction adopted by 
iGil dc Paz et alj 2003") , who derive their extinction from 
the Schlc gel et al, " ^998 ) pointing directly towards IC 10. 
That pointing may be confused by FIR emission from the 
galaxy itself. 

We adjust the Ha fluxes from IC 10 to account for 
1.1 magnitudes of internal extinction but this is a large 
source of uncertainty. The reddening extinction towards 
unem bedded stars is due ove rwhelmingly to Galactic 
dust llHunteHl200l Howev er. IVang fc Skillmal iTmgl 
and iBorissova et alJ l)2000() compare radio continuum 
and Br7 fluxes to Ha and find additional reddenings of 
E{B-V) ^ 1 {total E{B-V) ^ 1.5-2) towards embed- 
ded HII regions, implying extinctions as high as Aji 2.6 
towards these sources. Most of the Ha emission from 
IC 10 comes from regions that are at least somewhat 
embedded (see Figure EJ: 80% comes from lines of sight 
with H I columns above 10 M© pc~^ (1.25 x 10^^ cm^^) 
and the mean H I column associated with a bit of Ha 
emission ( ^ /(i/a)iV(H I)/ ^ /(77a) ) is 15 Mq pc'^. 
In the Milky Way, this column of gas (assuming it all lies 
between the obser ver and the H a) implies a reddening of 
E{B - y) ~ 0.4 l|Bohlin et al.l ll9781 and a correspond- 
ing i?-band extinction of ~ 1.0 magnitudes. Of course, 
the H I column is unlikely to all lie in front of the Ha 
emission, but there is a contribution from dust associ- 
ated with molecular g well. We use the Ha map to 
test the applicability of a "Schmidt Law" to IC 10 and 
this value is similar to the 1.1 m agnitudes of internal ex- 
tinction assumed bv iKennicuttI |l998.'l . so we adopt that 
value (1.1) for consistency but note it to be uncertain by 
^ 1 magnitude. 

3.3.3. VLA High Resolution HI Map 

We use the high resolution VLA H I map of 
iWilcots fc MilleJ ((122SI), reduced using natural weight- 
ing, which has a resolution of 11" (44 pc). The VLA map 
contains mostly data from intermediate arrays (B and C, 
minimum baseline ~ 35 m, with just 10 minutes of D ar- 



ray data) and is therefore not sensitive to spatial scales 
> 15'(although power on scales as small as ~ 7.5'may be 
attenuated by ~ 40%). In this work, we are interested 
in the H I over the central disk (where the CO is). In 
this region the structure of the H I varies on ~ 1' scales 
(the holes and filaments in Figure^ and we do not ex- 
pect significant loss of power on these scales. However, 
the H I surrounding IC 10 is very extended and spatial 
filtering is important to large scale studies. For example, 
IWilcots & Miller ( 1998) find only 60% of the fiux recov- 
ered by Shostak & Skilhnan ( 1989} using the WSRT. In 
turn, jghostak & Skillman (1989) recover only half of the 
value found by single dish telescopes. The extended emis- 
sion being resolved out by the interferometers has low 
column density (- 1 x lO^" cm-^ lWilcots fc Mill er' 19981. 
a factor of ~ 10 lower than t he typical colum ns associ- 
ated with CO detection (see iShostak fc S killman 1989|, 
for a map of extended emission missed by the WSRT). 
The large fiux discrepancies come from the large spatial 
extent of the H I and should only represent a ~ 10% 
correction on the compact, high column structures asso- 
ciated with the CO. 

3.3.4. K-hand Photometry 

IC 10 is part of the 2MASS Large Galaxy Atlas 
llJarrett et al. 2003). We use the iiT-band image to trace 
the stellar population in IC 10. Because IC 10 is at a low 
Galactic latitude (6 = —3.3°) foreground stars represent 
a serious source of contamination. In order to remove 
these foreground stars, we mask out the brightest pix- 
els in the 2MASS image. Our threshold for identifying 
"bright pixels" corresponds to a stellar surface density of 
1800 M© pc~^. We removed all pixels above this value 
from the image, as well as all data adjacent to these pix- 
els with a surface density of ~ 600 M© pc~^ (to ensure 
that we clipped the tail of the point spread function for 
the stars we remove). The highest stellar surface den- 
sity we find for the disk of IC 10 is < 500 M©, so the 
galaxy should be unaffected by the masking. After re- 
moving bright pixels, we apply a 15" median filter to 
the whole image, using the median values to replace the 
data we removed. Finally, we adopt a if -band mass to 
light rati o of 0.5 M© / L©.y , consistent with the results 
found by iSimon et al.l l)2005D in their study of rotation 
curves of dwarf galaxies. The scatter they find in their 
mass-to-light ratios is ~ 50% (ranging from 0.3 to 0.7 
Mq / I^q^k) so the stellar surface density is uncertain by 
the same amount. The resulting K-hand surface density 
map has a resolution comparable to the BIMA survey 
(because of the filtering) and is largely uncontaminated 
by bright foreground stars. Figure ^ shows this map as 
contours plotted over the extended H I distribution. 

3.4. Overview of GMC Property Measurements 

This section summarizes the algorithm described in de- 
tail by Rosolowskv & Lcrov (2005), which we use to mea- 
sure resolved GMC properties and correct these measure- 
ments for biases introduced by our limited resolution and 
sensitivity. We apply this algorithm to the BIM A survey 
and to the high resolution OVRO data sets from lWalter I 
( 2003) to produce the cloud property measurements in 
Tables H and 

First, we construct a mask containing all high signifi- 
cance signal in the data cubes. For the BIMA survey we 
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use the mask generated as described in ii3.1.2l and sim- 
ply assign emission to the nearest local maximum (see 
the cloud assignments in Figure^)). For the OVRO data, 
we include all regions with two adjacent velocity chan- 
nels both containing emission above 4(t intensity. We ex- 
pand the mask to include all emission with two adjacent 
velocity channels above 2a that is contiguous with the 
4(T peaks. We then identify significant, independent lo- 
cal maxima within each cloud. Here "significant" means 
that the maxima are at a significantly higher intensity 
(2(7 greater) than either the edge of the cloud or the high- 
est isosurface shared with other local maxima. An "inde- 
pendent" maximum is separated from all other maxima 
by at least a velocity channel or a full beam width. 

From the emission uniquely associated with each max- 
imum (i.e. within the lowest isosurface containing only 
that maximum), we measure the size, line width, and lu- 
minosity for that cloud. We make the measurements us- 
ing intensity-weighted moment methods (i.e. we measure 
spatial and velocity dispersions). From the measured dis- 
persions, w e calculate the r adius of the cloud using the 
definition of lSolomon et alJ |1987). The line width is the 
full-width at half-maximum of the integrated spectrum 
of the cloud. The luminosity is the integrated emission 
from the cloud. We correct these measurements for bi- 
ases due to the finite sensitivity and resolution of the 
astronomical data. We correct for the finite sensitivity 
by extrapolating the measured properties to those we 
would expect for a data set with perfect sensitivity (by 
fitting each property as a function of boundary isosurface 
value and extrapolating to a boundary of K). We cor- 
rect for the effects of beam convolution on the measured 
size of the GMCs by deconvolving the beam size from 
the measured size in quadrature (separately for the ma- 
jor and minor axes). Since the detailed description of the 
algorithm and its characterization i s beyond the scope of 
this p aper, we refer the reader to iRosolowskv fc Lerovl 
ll20(B . 

4. RESULTS 

4.1. Comparison of the Three CO Datasets 

How do the emission properties from the BIMA sur- 
vey, the ARO 12m data, and the OVRO observations 
compare? We convolved the BIMA survey, OVRO data, 
and the H I data (first clipped at Scr) to the resolution of 
the ARO 12m data. Figure shows these data for each 
ARO 12m pointing with detected emission (the H I spec- 
tra are arbitrarily normalized). Based on these spectra 
and Figure we draw the following conclusions: 

1. The agreement between the central velocities among 
all three datasets is excellent. Further, the velocities are 
consistent with the H I emission (dashed line in Figure 

Towards the brightest fines, the width of CO emission 
detected by the ARO 12m is comparable to, but always 
a bit smaller than, that of the H I emission. 

2. Both interferometers miss emission, and OVRO misses 
more than BIMA. This is probably a result of interferom- 
eters filtering out extended emission and not due to low 
signal-to-noise or to calibration errors. Both the BIMA 
survey and the OVRO data are more sensitive than the 
ARO 12m (both data sets have RMS noise < 10 mK 
over an ARO 12m beam at 6 km s^^ velocity resolution, 
compared to a typical noise of 15 niK for the ARO 12m 



spectra). Further, the OVRO data are missing more fiux 
than the BIMA survey, despite having a better signal-to- 
noise. Finally, the difference is not only a gain offset, as 
one would expect from a calibration discrepancy. Rather, 
the line widths of the interferometer data (particularly 
the OVRO data) are smaller than those found by the 
single dish, implying that an extended component with 
a larger velocity width contributes to the single dish data 
but not the interferometer data. 

How much emission is resolved out by the interferome- 
ters? On average, BIMA recovers 50% of the integrated 
intensity found by the ARO 12m, and OVRO recovers 
30%. BIMA finds a line width that is on average 90% 
of that recovered by the ARO 12m, while OVRO line 
widths are, on average, 60% that of the ARO 12m for 
the same pointing. We simulated BIMA D array obser- 
vations of several gaussian sources at the declination of 
IC 10. The observations recover 95% of the fiux for a 
20" (FWHM) source, 60% for a 30" source, 30% for a 
40" source, and 14% for a 40" source. If the CO struc- 
tures in IC 10 are ~ 35" (160 pc) we might expect to lose 
half of the fiux. The actual scales may be slighly more 
compact since there are inefficiencies associated with sig- 
nal identification and non-ideal observing conditions that 
will affect diffuse emission or long baselines. 

3. We detect CO emission with the ARO 12m only at 
pointings which also show emission in the BIMA survey 
(see point 4). Therefore the diffuse emission appears to 
be associated with the larger GMCs. Further, the ARO 
12m pointings cover almost all of the CO emission found 
in the BIMA survey. Only ^ 1 x 10^ K km s~^ pc^ 
(BIMA clouds 7, 12, 13, 14, and 15) of the emission from 
the BIMA survey (about 10% of the total) lies in clouds 
not targeted by the ARO 12m pointings. 

4. Two pointings deserve specific commentary: num- 
bers 7 and 14. Pointing 7 is detected by the ARO 12m 
and Figure El shows that it is also detected in the BIMA 
survey. However, the emission is not strong enough to 
be included in our mask. Pointing 14 is detected in the 
BIMA survey, but the ARO 12m emission appears to be 
contaminated by emission in the off position. 

Our comparison of the three data sets suggests that 
~ 50% of the CO emission from IC 10 comes from a 
high velocity width, spatially extended component that 
is resolved out by both OVRO and BIMA. OVRO further 
resolves out another 20% of the CO emission seen by 
the ARO 12m. Similar res ults have been fo und in the 
Milkv Way, M 31 and M 33 (|Po1k et a,1. 1 988: .BlitTjll 9851: 
IRosolowskv et alJl2003j) . In those galaxies, too, diffuse 
gas, or an extended grouping of small molecular clouds 
indistinguishable from diffuse gas, may contribute a large 
portion of the CO emission along a line of sight. 

4.2. The Total Content of Molecular Gas and CMC 
Properties 

Figure^ and Table ^ present the results of the BIMA 
survey. The total luminosity from all 10 CO detections 
with the ARO 12m of 1.6 ± 0.3 x 10^ K km s"! pc^ and 
this is our formal lower limit to the total CO luminosity of 
IC 10. If we stack all of the ARO pointings and integrate, 
the luminosity rises to 2.0 x 10^ K km s~^ pc^ with the 
increase of 25% due emission not detected in individual 
pointings (a signature of a low-intensi ty diffuse c ompo- 
nent, perhaps similar to that found bv llsrae il lT997l) . The 
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TABLE 4 

Cloud Properties from the BIMA Survey 



Clniirl dl 

V_>lULiVi. -H- 








<5 




^ctr 


^ FW H M 


T .1 ITTI inoQif""\r 

i_J LlllllllUol L y 


Mr 




^ J ZUUU ) 




(J2000) 




(km s ) 


(km s ) 


^lU iV Kill h pC j 


ivi0j 








(3) 








V>) 




(7\ 


Bl 


Oh 19m 


58.6s 


59° 


18' 


40.6" 


-367.4 


7.7 


30. 


131. 


B2 


Oh 20m 


0.8s 


59° 


19' 


4.9" 


-366.3 


12.4 


38. 


164. 


B3 


Oh 20m 


1.3s 


59° 


16' 


42.3" 


-364.0 


16.1 


49. 


213. 


B4 


Oh 20m 


11.9s 


59° 


20' 


2.2" 


-343.6 


13.6 


51. 


221. 


B5 


Oh 20m 


17.2s 


59° 


18' 


43.3" 


-335.9 


17.4 


50. 


218. 


B6 


Oh 20m 


19.6s 


59° 


20' 


48.1" 


-337.8 


10.5 


75. 


326. 


B7 


Oh 20m 


21.2s 


59° 


20' 


8.1" 


-330.6 


11.4 


41. 


180. 


B8 


Oh 20m 


21.6s 


59° 


17' 


8.8" 


-339.7 


10.5 


102. 


443. 


B9 


Oh 20m 


22.9s 


59° 


21' 


18.3" 


-329.6 


7.5 


122. 


533. 


BIO 


Oh 20m 


27.6s 


59° 


17' 


26.2" 


-329.6 


15.3 


57. 


246. 


Bll 


Oh 20m 


27.7s 


59° 


16' 


59.4" 


-330.9 


14.9 


238. 


1036. 


B12 


Oh 20m 


29.8s 


59° 


19' 


34.0" 


-326.4 


4.8 


12. 


52. 


B13 


Oh 20m 


31.3s 


59° 


18' 


9.9" 


-324.5 


12.0 


25. 


110. 


B14 


Oh 20m 


32.6s 


59° 


17' 


21.6" 


-311.1 


15.0 


11. 


48. 


B15 


Oh 20m 


34.4s 


59° 


19' 


24.0" 


-317.7 


21.4 


17. 


76. 


B16 


Oh 20m 


48.1s 


59° 


18' 


58.1" 


-330.5 


13.9 


61. 


266. 



TABLE 5 
CMC Properties in IC 10 



Cloud 


a 






<5 




■"ctr 


Radius 


Vfwhm 


Luminosity 


loglO ^Lum 


log 








(J2000) 




(J2000) 




(km s~^) 


(pc) 




(km s"-'-) 


(10^ K km s-i pc2) 


(Mq) 




(Mq) 




(1) 


(2) 




(3) 






(4) 


(5) 




(6) 


(7) 


(8) 




(9) 


(10) 


Bl 


Oh 19m 


58.6s 


59° 


18' 


40.3" 


-366.6 


11.8 ± 


6.1 


4.0 ± 1.0 


13.4 ± 2.8 


4.8 ± 0.1 


4, 


,6 ± 0.2 


0.6 ± 0.5 


B2 


Oh 20m 


0.9s 


59° 


19' 


2.2" 


-367.4 


22.8 ± 


9.6 


5.7 ± 2.0 


9.6 ± 7.0 


4.6 ± 0.2 


5, 


1 ± 0.3 


3.3 ± 3.8 


B4 


Oh 20m 


12.0s 


59° 


20' 


2.5" 


-342.7 


29.4 ± 


7.2 


5.6 ± 1.0 


20.2 ± 2.7 


4.9 ± 0.1 


5, 


,2 ± 0.2 


2.0 ± 1.1 


B5 


Oh 20m 


17.3s 


59° 


18' 


42.0" 


-334.0 


26.2 ± 


5.5 


6.2 ± 1.3 


25.3 ± 4.7 


5.0 ± 0.1 


5, 


,3 ± 0.2 


1.7 ± 1.0 


B6 


Oh 20m 


19.6s 


59° 


20' 


48.0" 


-338.2 


17.7 ± 


3.9 


6.0 ± 1.2 


27.5 ± 5.3 


5.1 ± 0.1 


5, 


1 ± 0.2 


1.0 ± 0.5 


B8 


Oh 20m 


21.7s 


59° 


17' 


9.6" 


-339.5 


24.1 ± 


2.6 


6.1 ± 0.6 


69.5 ± 3.2 


5.5 ± 0.0 


5, 


,2 ± 0.1 


0.6 ± 0.1 


B7 


Oh 20m 


22.1s 


59° 


20' 


4.9" 


-329.3 


18.3 ± 


6.3 


4.0 ± 1.1 


18.0 ± 5.4 


4.9 ± 0.1 


4. 


,7 ± 0.2 


0.7 ± 0.5 


B9a 


Oh 20m 


22.3s 


59° 


21' 


5.6" 


-330.5 


15.9 ± 


4.6 


3.0 ± 0.8 


14.1 ± 3.2 


4.8 ± 0.1 


4. 


,4 ± 0.2 


0.4 ± 0.3 


B9b 


Oh 20m 


22.5s 


59° 


21' 


21.4" 


-328.7 


24.5 ± 


6.3 


4.3 ± 1.1 


18.7 ± 5.3 


4.9 ± 0.1 


4. 


,9 ± 0.2 


1.1 ± 0.7 


B9c 


Oh 20m 


23.7s 


59° 


21' 


17.9" 


-333.8 


27.2 ± 


7.7 


5.5 ± 2.5 


22.4 ± 5.8 


5.0 ± 0.1 


5, 


,2 ± 0.3 


1.6 ± 1.5 


Blla 


Oh 20m 


27.2s 


59° 


16' 


53.8" 


-334.0 


22.1 ± 


7.3 


10.4 ± 2.5 


76.4 ± 24.1 


5.5 ± 0.1 


5, 


,7 ± 0.2 


1.4 ± 1.0 


Bllb 


Oh 20m 


27.3s 


59° 


17' 


5.5" 


-331.1 


14.9 ± 


3.0 


7.4 ± 1.1 


52.5 ± 8.5 


5.4 ± 0.1 


5, 


,2 ± 0.1 


0.7 ± 0.3 


Bllc 


Oh 20m 


28.1s 


59° 


16' 


57.0" 


-325.1 


17.8 ± 


4.5 


11.2 ± 2.2 


64.3 ± 12.0 


5.4 ± 0.1 


5, 


,6 ± 0.2 


1.5 ± 0.8 


Blld 


Oh 20m 


29.0s 


59° 


17' 


4.6" 


-327.4 


18.4 ± 


3.9 


7.9 ± 1.5 


25.6 ± 4.8 


5.0 ± 0.1 


5, 


,3 ± 0.2 


1.9 ± 1.0 



ARO 12m pointings overlap ^ 90% of the emission found 
in the BIMA survey, so we estimate that the total CO lu- 
minosity of IC 10 is - 2 X 10^ X gig = 2.2 x 10^ K km s^^ 
pc^. This number is quite uncertain because our mask- 
ing algorithm is chosen to avoid false positives (rather 
than for completeness). We estimate an upper limit by 
noting that the inclusion of 11 pointings without indi- 
vidual CO detections raises the integrated luminosity by 
4 X 10^ K km s^^ pc^. The total area in IC 10 with an 
H I surface density > 10 Mq pc^ (and therefore likely to 
harbor molecular gas) corresponds to ^ 16 times the area 
of the 12m beam. Of this area about 6 times the area of 
the 12m beam is already covered by our ARO observa- 
tions. Four of our nondetection pointings are extremely 
unlikely sites for molecular gas emission. Therefore, we 
might expect another ^ x (4 x 10^) 6 x 10^ K km s~^ 
pc^ in diffuse emission. We therefore suggest 2.8 x 10^ K 
km s~^ pc^ as an upper limit to the CO luminosity. The 
BIMA survey recovers about half of our best guess at the 
luminosity — the 16 GMCs listed in Table0|have a total 



luminosity of 1.0 x 10^ K km s^^ pc^ (which includes a 
^ 30% sensitivity correction as described in t i3.4l) . 

We used the algorithm of Rosolowsky & Leroy (;2005f l 
summarized in t|3.4l to identify clouds an d then measur e 
their properties using the OVRO dataset ijWalter 12003(1 . 
Table gives our measurements of the properties of 14 
GMCs. Column (1) gives the cloud name (which also 
indicates the BIMA survey cloud to which the OVRO 
cloud most nearly corresponds); columns (2) and (3) give 
the intensity-weighted average position of the cloud; col- 
umn (4) gives the intensity-weighted average velocity of 
emission from the cloud; column (5) gives the equivalent 
spherical radius of the cloud ( Solomon et al. 1987) in par- 
sees, after correction for the effects of beam convolution; 
column (6) gives the FWHM line width of the cloud; col- 
umn (7) gives the cloud luminosity in K km s~^ pc^; 
column (8) gives the mass derived from this luminosity 
assuming a Galactic CO-to-Il2 conversion factor; column 
(9) gives the dynamical mass of the cloud, calculated as- 
suming virial equilibrium and a p oc density profile. 
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Fig. 6. — ARO 12m, BIMA, and OVRO spectra compared for common pointings which show emission. The BIMA and OVRO data 
have been convolved to the resolution of the ARO 12m. The H I spectrum, normalized to 80% of the y-axis range, is overplotted as a 
dashed line to indicated the velocity width of the H I. 



SO that Myir = 189 R Vp-^jjj^.j, with R and Vfwhm 
as defined above; column (10) gives the ratio of virial 
to luminous mass for the cloud. The virial masses are 
comparable to the luminous masses, suggesting that the 
CO-to-H2 conversion factor does not vary significantly 
from the Galactic value. Figure [7| shows a map each 
of cloud. Ellipses indicate the size of each cloud before 
any corrections for resolution or sensitivity effects are ap- 
plied, so that the sizes shown is directly comparable to 
the structures in the data. In Table[5land in the analysis 
below, however, we use the corrected values. 



We include a brief word of warning regarding cloud 
names. For clarity we refer to clouds measured in the 
OVRO data using the same system we used for the BIMA 
survey, using "a,b,c" when a BIMA cloud is resolved into 
several GMCs by OVRO. However, we have not tried to 
systematically identify a single set of clouds using both 
data sets, and the association implied by a shared name 
should not be weighted too heavily. For example, though 
no cloud BIO exists in the OVRO data, emission from 
that object certainly corresponds to some of the emission 
found in clouds Blla,b,c, and d. 
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4.3. Comparison to CMC Properties From the 
Literature 

There are several observations of IC 10 GMCs in the 
literature: do the properties we measure agree with these 
studies? Table El shows a comparison of the properties of 
three GMCs in the bright CO complex in the southeast 
part of the galaxy (cloud 15 in Figurc01 ARC 12m point- 
ing 6 in FigureEJ . TableElshows the prop erties measured 
by two previous studies (fWilsonl995: Ohta et alJll992(l 
and this work. All three works decompose the emission 
in essentially the same manner, though some of the vari- 
ations in Table may arise from decisions about which 
emission to assign to which CMC. To properly compare 
these observations we have adjusted the sizes measured 
by the two previous studies to match our definition of the 
radius, our adopted distance, and our adopted Galactic 
CO-t o-H2 convers ion factor — both Ohta et al. ( 199^ 
and Wilson ( 1995) quote full width half maximum sizes 
and th ey adopt di s tances of 1.3 and 0.82 Mpc, respec- 
tively. lOhta et all ljl992|) labels the GMCs in q uestion 
'NCI,' 'NC2,' and 'SC (1, 2, and 3 in Table Wilsoul 
(p995) calls them 'MCI,' 'MC2,' and 'MC3;' they are 
'Bllb,' 'Bllc,' and 'Blla,' respectively, in our Table [3 
The sizes agree well for two of the three clouds, and we 
measure a notably lower size for CMC 1; there is a sig- 
nificant (~ 50%) scatter i n the line widths. The fluxes 
measured bv lWilsonI l)1995f) are systematically lo wer than 
what we measure — integrated over the complex. iWilsoiil 
lfT99,5.) measures half of our flux. The discrepancy in 
fluxes and ou r 15% larger adopted distance explain why 
iWilsonI l|1995|) derives a higher CO-to-H2 conversion fac- 
tor for the GMCs than we do. The combination of a 
larger adopted distance and the systematically higher 
fluxes found by Ohta et al. ( 1992) explain the similar- 
ity in the derived CO-to-H2 conversion factors. 

Which measurements are closer to the true fluxes? The 
complex in question corresponds to pointing 6 with the 
ARO 12m (see Figures IH and and Table 0) and cloud 
BIMA 11 in Table H BIMA and OVRO find similar 
fluxes for the region (both corresponding to « 10® M© for 
our adopted Xco) and the ARO 12ni recovers about 2.5 
times this luminosity. IWilsonI l)1995t) finds fluxes lower 
than those found by BIMA, the newer OVRO dataset, 
or the ARO 12m. The NMA dataset appears consistent 
with the BIMA and OVRO result s (also fi nding « 10® 
Mq for our Xco)- Therefore the IWilsonI l|1995li fluxes 
represent an outlier from the other interferometric data. 
All four interferometric data sets resolve out a large frac- 
tion of the flux (as measured by the ARO 12m). 

5. DISCUSSION 

5.1. Molecular Gas Fraction and Depletion Time 

Using a Galactic CO-to-H2 conversion factor, 2 x 10^" 
cm~^ (K km s~^)~^ (appropriate for the GMCs but per- 
haps not for diffuse gas, see t|5.2|l . the luminosity we es- 
timate for IC 10 translates to a molecular gas mass of 
9 X 10® Mq. This mass is small compared to the other 
components of IC 10 — the stellar, H I and dynami- 
cal masses are ~ 4, 2, and 15 x 10® Mq respectively 
— and this implies that only 4% of the gas is molecu- 
lar. For a star formation rate of ~ 0.2 Mq yr~^, the 
depletion time for the molecular gas associated with CO 
in IC 10 is 4 X 10^ years, compared to the median 



depletion time of 2 x 10^ years found in nearby dwarf 
galaxies (|Lerov et al.l 12005'!. The molecular gas mass 
is only about ~ 2% of the stellar mass. This value 
tends to be larger ^ 15%, in LMC-type d warf galaxies 
l|Young fc S covillc 1 99HlLerov et aH'mi^. In fact, the 
amount of molecular gas per stellar luminosity tends to 
be fairly constant among all star forming galaxies, large 
and small (in the i^T-band this ratio is ~ 0.07 Mq / 1uk,q 
iLerov et al.ll2005D . Although IC 10 is CO-deficient com- 
pared to larger galaxies, i t is CO-rich when compared to 
the SMC or NGC 1569. iMizuno et al. (2001) find that 
the CO luminosity of the SMC is « 8 x 10"^ K km 
pc^, more than an order of magnitude fainter than the 
CO luminosity of IC 10, while NGC 1569, which has a 
CO luminosi ty of ~ 10^ K km pc^ ijCreve et al.ll996t 
iTavlor et al.]ll999l) . This yields an even shorter molecu- 
lar gas depletion time for these two systems than what 
we observe in IC 10. 

5.2. The CO-to-H2 Conversion Factor in IC 10 

With lower abundances of carbon and oxygen, harder 
radiation fields, and less dust to shield each parcel of 
gas, dwarf galaxies may be expected to display a dif- 
ferent relationship between CO emission and molecular 
hydrogen content. Calibrating the CO-to-Il2 conversion 
factor, Xcoi has often been a goal of CO studies of dwarf 
galaxies fe.g. IWilsonill99,'TD . The key to measuring Xco 
in other galaxies is to find an independent method of 
measuring the amount of molecular gas present. We at- 
tempt an independent estimate of the mass by measuring 
the size and line width of a molecular cloud from its CO 
emission and then calculating the dynamical mass under 
the assumption of virial equilibrium. 

In our analysis of the high resolution CO data, we find 
that the IC 10 clouds are indistinguishable from GMCs 
analyzed in the same way in M 31 and M 33 and that they 
are very similar to Milky Way clouds. If the IC 10 clouds 
are virialized, then the mean C0-to-Il2 conversion fac- 
tor in the CO peaks is 2.6 x 10^° cm-^ (K km s''^)-'^ 
(if the clouds are only marginally bound then Xco 
will be half of this value). Virial mass studies in the 
Milky Way yield a CO-to-H a conversion fac tor Xco ~ 
3 X lO^o cm-2 (K km s~^)~^ jSolomon et alJlT987|) . sim- 
ilar to the one we measure in IC 10 within the uncer- 
tainties; interferometric studie s of M 31 an d M 33 yield 
app roxim ately the sa i ne res ult iRosolowskv et al. ( 2003); 
iRosolowskv fc Lerovl l|2005j) . Because, in the Milky Way, 
the Xco value derived from gamma-ray observations is 
thought by many to be more reliable (it is indepen- 
dent of the dynamical state of the GMCs), we adopt 
the Galactic value of Xco ob tained by those studies, 
2 X 10^° cm-^ (K km s-^"^ llStroug fc Mattoxl FmM 
iDame et al.l l200lf) in the remainder of this study. 

In igllwe found that OVRO may resolve out 70% of 
the emission. The measurements of the CMC properties 
are probably robust, since the GMCs are compact rela- 
tive to the ~ 100 pc scales we expect OVRO to resolve 
out and the OVRO data have good sensitivity. How- 
ever, the CMC properties measured from the OVRO data 
do not constrain the CO-to-Il2 conversion factor in the 
extended emission. One possibility is that the CO re- 
solved out by OVRO and BIMA comes from a spatially 
extended collection of small GMCs. Below we find evi- 
dence for a CMC mass spectrum with a power law index 
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Fig. 7. — Intensity maps of GMCs in IC 10. The figure shows all CO emission near the center of the GMC. The circles indicate the 
measured sizes of the GMC sensitivity or resolution corrections. At the distance of IC 10 (950 kpc), 1" corresponds to 4.6 pc, so the 40" 
boxes shown have spatial sizes of 92 pc. 



of —2.0, which implies that there may be as much mass in 
GMCs below our completeness limit as above it. If these 
low mass GMCs make up the extended structure that is 
resolved out by BIMA and OVRO, then we expect that 
the CO-to-H2 conversion factor from the OVRO clouds 
will apply, at least approximately, to all of the CO emis- 
sion. 

B ased on observations of the 158 /im [CII] emission 
line. lMadden et aP l)1997!) suggested a CO-to-H2 conver- 
sion factor much higher than the Galactic value. They 



mapped IC 10 at ~ 1' resolution and found that it is 
luminous in the [CII] line compared to the CO luminos- 
ity. In the northern and western regions, they found that 
the minimum amount of hydrogen needed to produce the 
observed [CII] emission implies a substantial mass of gas 
that was not inferred from the H I or CO emission. They 
suggested that in parts of IC 10 molecular hydrogen may 
exist under conditions of low extinction. Ay ~ 1 ~ 2, 
so that CO is dissociated but self-shielded II2 is abun- 
dant. They argue that in these regions the II2 column 
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TABLE 6 

IC 10 SE GMC Property Comparison 



GMC 


Property 


Ohta et al. (1992) 


Wilson (1995) 


This Paper 


1 
± 


ivaclius (pc) 


OQa.b.c 


Qoa.b.c 


10 


1 


Vfwhm (km s-i) 


7.3 


10.9 


7.4 


1 


Luminous Mass (10^ Mq) 


290''''' 


190"'" 


250 


2 


Radius (pc) 


iga.b.c 


23a,b,c,d 


18 


2 


Vfwhm (km s'l) 


7.3 


8.9 


11.2 


2 


Luminous Mass (10^ Mq) 


230"'" 


120"'" 


250 


3 


Radius (pc) 


24a, b,c 


24a,b,c,d 


22 


3 


Vfwhm (km s'^) 


< 8.2 


12.8 


10.4 


3 


Luminous Mass (10^ M0) 


490"'" 


130"'" 


320 



"Value has been adjusted to our assumed distance of 950 kpc. 
Geometric mean of 2D size from the literature. 
"Adjusted to our radius definition from FWIfM. 
'^ Cloud unresolved in one dimension so size is a maximum. 
"Adjusted to our adopted value of Xco. 



may exceed the H I column by a factor of 5. Although we 
a dopt a Galactic Xco, we can not rule out the evidence 
of IMadden et all l|1997j) for a large reservoir of non-CO- 
emitting molecular gas outside the central GMCs. We do 
not account for such gas in our discussion of star forma- 
tion efficiencies because such gas must be warm, diffuse, 
poorly shielded, and, as such, does not seem to be a likely 
locale for star formation (the excitation temperature of 
the 158 /im [CII] li ne is 92 K and th e CO-free molecular 
gas posited by Mad den et alJ l)1997D exists at extinctions 
of Ay ~ 1 - 2 magnitudes). Further, if such gas ex- 
ists in other galaxies it will be similarly unaccounted for 
by the CO luminosity and should thus be left out for a 
self-consistent comparison. 

There may be evidence for a reservoir of warm, CO- 
free H2 beyond the GMCs, however observations do not 
suggest the presence of a hidden reservoir of cold gas. In 
a stud y of the dust continuum in the southeast part of 
IC 10. iBolatto et all ((20000 considered and rejected the 
possibility of a large reservoir of cold molecular gas in 
that region. Although they find an excess of long wave- 
length infrared emission, they cite the lack of CO self 
absorption and the normal CO (2^1) to CO (1^0) 
ratios as evidence that the long wavelength emission is 
not indicative of a massiv e reservoir of cold molecular 
gas. iThronso n et al.l (|199(J') also finds the amount of 155 
/im dust emission in IC 10 to be consistent with the H I 
emission from the galaxy. A large population of hidden 
molecular gas is not necessary to explain their observa- 
tions, though IC 10 does show a mild excess at 155 /^m. 

For the rest of this paper, we adopt a CO-to-H2 con- 
version factor of 2 X 10^" cm^^ (K km s^^)^^ for all of 
the CO emission. We are confident in the applicability of 
this value of Xco to the GMCs but less certain whether 
it is appropriate for diffuse CO emission. We neglect 
the possibility of warm H2 untraced by CO because we 
have no observational handle on such gas and it seems 
unlikely to form stars, but we emphasize that any such 
excess of molecular gas must exist outside the most mas- 
sive GMCs or we would see evidence for it in the virial 
masses we measure from the OVRO data. We note the 



following conversions for our adopted Xco ■ an integrated 
intensity of 1 K km s^^ corresponds to a molecular gas 
surface density of 4.4 Mq pc~^, including helium, and 
therefore a luminosity of 1 K km s~^ pc^ translates into 
a molecular mass of 4.4 Mq. 

5.3. Mass- Size- Line Width Relations in IC 10 

In the Milky Way, M 33, and M 31 molecular gas is 
concentrated in GMCs which exhibit s c aling laws relating 
their properties iSolon ion et al.lll987t iRosolowsk v et all 
120031 iRosolowskv 200^. Thes e relat ionships, often re- 
ferred to as "Larson's Laws," l)Larsonlll9 81') relate the 
size of a GMC to its line width, the CO luminosity 
to the dynamical mass, and the luminosity to the line 
width. In this section we compare clouds in IC 10 to 
clouds from the Milky Way, M 33, and M 31. M 31 
and M 33 are at distances comparable to that of IC 10, 
< 1 Mpc, and the interferometric data used for the com- 
parison have similar spatial resolutions (« 20 pc). We 
measure their properties using the same algorithm used 
to analyze the IC 10 clouds (a reanalysis of the dat a 
presented in IRosolowskv et ani2003H Rosolow skv 112005)) . 
Thus the M 31 and M 33 data should represent an excel- 
lent "control" sample, and any differences between GMC 
properties in these systems and IC 10 should be a re- 
sult mainly of environmental factors, not observational 
or analytical biases. The Milky Way data have not been 
analyzed in the same manner as the other data sets — 
they just cons i st of the GMC properties measured by 
iSolomon et alJ l|1987j) — so systematic differences may 
bias the comparison. However, we have attempted to 
correct for sensitivity and re solution biases in our data 
and the ISolomon et aP l)1987(l data have good sensitivity 
and spatial resolution (being Galactic data) so we antic- 
ipate the magnitude of these biases to be small. 

Figure IHl shows that clouds in IC 10 lie on or near the 
scaling laws for GMCs in the Milky Way, M 31, and 
M 33. This conclusion is reinforced by Figure O which 
shows three physical parameters — the surface density, 
the virial to luminous mass ratio, and the scatter about 
the size-velocity dispersion relation {dy cx R'^-^) — for 
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clouds in IC 10, the Milky Way, M 31, and M 33. Table 
dshow the results of two sided KS-tests comparing these 
physical parameters between IC 10 and the Milky Way 
(column 2) and M 31+M 33 (column 3). The distribu- 
tion of physical parameters for IC 10 GMCs agrees very 
well with the distribution found for the M 31 and M 33 
GMCs (< 2(7 difference for each parameter). Since the 
M31/M33 CMC properties are also derived from inter- 
ferometric data and are the result of the same analysis 
used to produce the IC 10 CMC properties, we attach 
particular weight to this comparison. The KS test de- 
tects ~ 3(7 differences between in IC 10 the Milky Way in 
the both surface density of clouds and their scatter about 
the size- velocity dispersion relation. Figure shows that 
these differences are relatively small however and since 
the same differences exist between the Milky Way data 
and the M 31/M 33 data they may represent a system- 
atic difference stemming from either observational biases 
or differences in the measurement algorithm. The CO- 
emitting clouds resolved by OVRO in IC 10 are very sim- 
ilar to GMCs in the Milky Way. They are indistinguish- 
able from GMCs in M 31 and M 33 observed and analyzed 
in the same manner as the IC 10 clouds. 

5.4. The Mass Spectrum of GMCs m IC 10 

Do GMCs in IC 10 also exhibit the same distribution 
of masses as GMCs in other galaxies? We calculate the 
cumulative distribution function (CDF) for CMC masses 
in our data (calculated from the CMC luminosities), so 
that the value of the CDF at a particular mass is equal to 
the fraction of clouds with masses greater than or equal 
to that mass. The power law index of the CDF indicates 
how "top heavy" or "bottom heavy" the distribution of 
cloud masses is. Figure ^| shows the CDF for clouds 
from the BIMA survey and GMCs above 5 x 10^ Mg 
from [Solomon ct al. (1987). The best fit power law to 
the BIMA survey CDF yields ^ cx M-2 0±0-2, a spec- 
trum with approximately equal mass in each logarithmic 
bin. This index agre es with the Milky Way slope we find 
when we consider the lSolomon et aD l)1987j) data over the 
same range of masses (— 1.85±0.05), but is steeper (more 
"bottom heavy" ) than the slop e of w —1.6 derived from 
all of the iSolomon etaP l)1987() data. The OVRO data 
also yield a power law index of ~ —2 (though they do not 
represent a complete sample and show a normalization 
offset). The exact power law index of the mass distribu- 
tion in IC 10 is quite uncertain because it depends heavily 
on the identification of GMCs in the BIMA survey (for 
example the large cloud in the southeast is resolved into 
three separate, smaller clouds, by OVRO). By altering 
our decomposition of the survey, we are able to obtain 
power law indices between — 1.8±0.2 and — 2.2±0.2. The 
gray region in Figure 1101 shows uncertainties associated 
with counting errors and uncertainties in the mass, but 
not decomposition. 

We have tested for consistency by realizing 1000 test 
samples of Milky Way GMCs and comparing their mass 
distribution to that of the IC 10 GMCs. Each test 
sample contains 16 clouds (to match the IC 10 sample) 
rando mly drawn from the population of ISolomon et ahl 
(jl987|) clouds with masses > 5 x 10'* M©, allowing re- 
peats. We compare each test sample to the population 
of clouds in IC 10 using a two sided KS test. As a control, 
we perform the same test using pairs of test samples from 



the Milky Way data. The median comparison of IC 10 
GMCs to Milky Way GMCs (over the 1000 test samples) 
showed more difference than 65% of the control compar- 
isons. Any differences between the two populations are 
thus of only ^ la significance. 

A mass distribution with a power law index of —2 im- 
plies a significant amount of gas below our completeness 
limit. As mentioned in the previous section, these small 
clouds might explain the discrepancy between the single 
dish and interferometer results if they exist in an ex- 
tended layer (say, throughout the H I filaments) near the 
large GMCs. 

5.5. CO and HI 

The GMCs in IC 10 are found almost exclusively in 
high column density atomic gas filaments. Figure 0] 
shows that we find molecular gas only where we find 
atomic gas, usually above a surface density of T^hi = 10 
Mq pc-2 (iV(HI) = 1.25 X 10^1 cm-2). Figure ITTl shows 
the relationship quantitatively, plotting the fraction of 
the total CO emission as a function H I column density 
along the line of sight (in black) and the fraction of lines 
of sight with the specified column density that have asso- 
ciated CO emission (in gray). Half of the molecular gas 
emission comes from regions with H I column densities 
above 15 Mq pc~^ and 85% is found above a contour of 
'Shi — 10 Mq pc^^. A high column of atomic gas is not 
a sufficient condition, though, as only 30% of the area in 
IC 10 with > 10 Mq pc~^ has associated molecu- 
lar gas and only half of the lines of sight with S^f/ — 16 
Mq pc-2 {N{R I) = 2.0 X 10^1 cm-2) have associated CO 
emission. This discrepancy (necessary, but not sufficient) 
arises mostly as a result of the large region of relatively 
high column atomic gas to the east of IC 10. This high 
column gas is mostly dev oid of CO emission, har boring 
only one molecular cloud. lEneargiola et alJ lj2003|) found 
a similar result in their survey of M 33 — the CO emis- 
sion comes almost exclusively from within the H I fila- 
ments but the presence of a filament does not necessarily 
imply the presence of CO emission. Similar results can 
be seen com paring CO emiss ion from the LMC to H I 
l)Fukui et al.lll999t iKim et al.lll998() . Broadly, this is the 
same effect seem in many disk galaxies: the H I extends 
far out into the disk while molecular gas a nd star forma- 
tion are relatively centrally confined (e.g. iWong fc Blitd 
l20f)l. 

IWong fc Bm2 l|200^ and IBlitz fc Rosolowskvl pOOl 
have suggested that the hydrostatic gas pressure, Ph, 
can predict the ratio of molecular to atomic gas, f„ioi — 
Y^moi/'^Hi, over a region of a galaxy. The hydrostatic 
gas pressure may trace the volume density of gas, Pgas, 
since Ph = pgas Vg and the gas velocity dispersion, Vg, is 
often quite constant l)Blitz fc Rosolowskvl >2004". and ref- 
erences therein). The volume density of gas should be 
more relevant to the formation of II2 from H I than the 
surface density. We test whether these arguments hold 
in IC 10 by calculat i ng Ph using the formula derived by 
IBlitz fc Rosolowskvl l|2004j) for a stellar-dominated disk, 

P. = 0.84(GS.)°-^S,^£2_ (4) 

where Eg is the total surface density of the g 

the surface density of stars, Vg is the velocity dispersion 
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TABLE 7 

KS Test Comparisons of IC 10 GMCs 



Property 


vs. Milky Way 


vs. M 31+M 33 


Surface Density 


0.02=' 


0.65 


Size "'^ / Velocity Dispersion 


1 X 10-3 


0.41 


Virial Parameter 


0.25 


0.67 



^Entries are the probability of generating the observed dif- 
ference or greater from the same parent population (i.e. ran- 
domly) . 
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Fig. 8. — Larson's laws for clouds in IC 10 (black circles), compared to the Milky Way (black dots) and the Local Group spirals M 31 
and M 33 (gray circles), (a) The size- velocity dispersion relation, (b) The virial mass-luminosity relationship, (c) The luminosity-line 
width relationship. All four galaxies display very similar behavior, obeying the same scaling of their global GMC properties. 



of the gas, and ft„ is the seale height of the stellar disk. 
From the H I cube (Wilcots & Miller 1.99S,) , we measured 
the median Vg across the disk (for lines of sight with non- 
negligible H I emission) to be « 7.5 km s""'^. We assume 
the scale height of stars t o be comparable to that o f disk 
galaxies (~ 300 pc, see iBlitz fc R.osolowskvl l2(10^ and 
calculate fmoi using Xco = 2x 10^° cm~^ (K km s~^)~^ 
(appropriate to the GMCs but perhaps not diffuse gas, 
see 



Figure ^1 shows reasonable agreement between 
the IC 10 d ata and an extrapolation from the 
IBlitz fc R.osolow skv (2005) galaxies to low pressure. Ph 
and fmoi have a rank correlation coefhcient of 0.7 ± 0.1. 
The median ratio, Pq = is 7 x 10"* ks and is un- 

certain by a factor of two. This value of Pq (the pres- 
su re for wh ich f„ioi = 1) is consistent with the results 
of lBlitz fc R osolowskv (2005), who find Pq = 4 x 10"* ks 
and a nearly linear best fit relation. On average IC 10 is 
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Virial Parameter (log,o U„yUt^„) 

Fig. 9. — Distributions around the Larson relationships for clouds in IC 10 (hashed), the Milky Way (thick line), and M 31 and M 33 
(gray histogram), showing the similarities among the GMC populations in the four galaxies. 



poorer in molecular gas than one might expect by about 
a factor of two (it lies ju st under the best fit line from 
iBhtz fc Rosolowskvl2005|) . Since we used the BIMA sur- 
vey in this comparison, the discrepancy might be almost 
completely negated by including the resolved out flux. 

Table [S] shows that the pressure is a better predictor 
of fmoi than either the atomic gas surface density or the 
stellar surface density. The rank correlation between Ph 
and fmoi is 0.7 ± 0.1, higher than the rank correlation 
between fmoi and either "Ehi or 12^:. Both T.hi and 
are very highly correla ted with Smn/, which demon strates 
that the conclusion of iBlitz fc Rosolowskvl l)2005(l seems 
to hold in IC 10: H I is a necessary prerequisite for the 
presence of molecular gas, but not the best predictor of 
the ratio of molecular to atomic gas. The hydrostatic 
gas pressure offers an improved prediction of the ratio of 
molecular to atomic gas because the gravitational influ- 
ence of a stellar potential well is necessary to enable the 
formation of H2 out of H I filaments. 

5.6. Star Formation and Gas in IC 10 

We showed above that GMCs in IC 10 are similar to 
those in the Milky Way, M 31, and M 33. Does the molec- 



ular gas also form stars at the same rate as molecular gas 
in spirals? In this section we compare the star formation 
efficiency (SFR per unit gas) in IC 10 to that in larger 
galaxies. We perform these comparisons using the H I, 
CO, and Ha maps convolved to a common spatial resolu- 
tion of 275 pc (60", with 30" pixels). This is comparable 
to the scale height of the stellar disk in a spiral galaxy, so 
we expect the averaging to occur over roughly the same 
distance along the line of sight and perpendicular to the 
line of sight. We limit ourselves to regions of the stellar 
disk of IC 10 with stellar surface densities in excess of 
~ 50 M0 pc~^ (obtained from the if-band light), be- 
cause outside this region the three dimensional structure 
of the H I envelope is uncertain and it is unclear how to 
interpret the line of sight surface densities. 

For galaxies the size of the LMC or bigger, the 
efficiency with which stars form from molecular gas 
(and its inverse, the molecular gas depletion time) 
depends wealky on galax y mass and Hubble type 
llYoung fc Scovillefl 991: You ng et al.ll995HMurgia et alJ 
2002; Lerov et al. 2005). Figure [T5I shows the relation- 
ship between surface density of star formation and molec- 
ular gas surface density for a range of galaxies. Both 
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TABLE 8 

Rank Correlations with Molecular to Gas Content 



Property Rank Correlation with S^oj Rank Correlation with Ti^^l/'^HI 

T,Hi 0.6 ±0.1 0.4 ±0.1 

S, 0.4 ±0.1 0.5 ±0.1 

Ph 0.8 ±0.1 0.7 ±0.1 
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Fig. 10. — The cumulative distribution function for GMCs in the 
Milky Way (solid line) and IC 10 (black circles). We include only 
Milky Way clouds with M > Mq (the approximate sensitivity 
of our survey). The gray region shows the area of Icr uncertainty, 
but systematic effects from the identification of clouds are also 
more important and are not shown. 



Fig. 12. — The molecular-to-atomic ratio plotted as a function 
of midplane pressure in IC 10. Results and the best fit relation 
for spiral galaxies from Blitz & Rosolowsky (2005) are shown for 
comparison. We bin data with similar pressures along lines of sight 
with S* > 50 M0 pc~^, where the disk approximation may be 
valid, and apply an inclination correction of 0.67. 
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Fig. 11. — The fraction of CO emission from IC 10 above a given 
HI column density, is shown in black. The fraction of lines of 

sight with a given HI column density that also show CO emission 
is shown in gray. Most of the emission comes from regions of the 
galaxy with column densities above 10 Mq pc~^, and 50% of the 
emission arises from regions with ^ 16 Mq pc~^. Only 50% 

of lines of sight with this column density has associated molecular 
gas, however. Although T,hi ^ 10 is a necessary condition to find 
molecular gas, it is not sufficient on its own. 



dwarf galaxies and spirals obey a power law relationship 
of roughly 

^SFR = lU l^Mol ' (5) 

where T,sfr is the star formation surface density in units 
of Mq yr~^ kpc~^ and T,moI is th e molecular gas sur- 
face density in units of Mq DC~^ l)Murgia et ^ 120021 
Eerov ct al. 2005). Figure E| also shows that IC 10 very 
clearly does not fall on this trend. Rather, the data for 
IC 10 show much larger rates of star formation per unit 
molecular gas than is found in other galaxies — a median 
factor of 270 higher. IC 10 is not tmique in this regard: 
both the SMC (Mizuno ct al. 2001: Wilkc ct al. 200|) 
and the nearby starburst NGC 1569 also show higher 
SFR surface densities than their CO content would sug- 
gest unless Xco is considerably larger than Galactic in 
these sources. 

The total gas surface density, rather than the II2 sur- 
face den sity, is often used as a predictor of the star for- 
mation. iKennicutd l|1998D found that a single power law 
described the relationship total gas surface density and 
the star formation surface density across a wide range of 
galaxies. Data from that paper are plotted along with gas 
surface densities and Ha derived star formation surface 
densities (SFSDs) in Figure Figure ITU shows that in 
IC 10, the highest gas surface densities do roughly corre- 
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spond to the highest star formation rates, although the 
agreement with Kcnnicutt (1998) is poor. In particular, 
regions of high SFSD seem to have lower gas surface den- 
siti es in IC 10 than galaxies with the same SFSD from 
the Kennicutt (1998) sample. Furthermore, the scatter 
in SFSD for a given total gas surface density is very large, 
particularly around J^hi+H2 ~ 7 Mq pc~^. 

How can we explain the extremely high efficiency of 
star formation in IC 10 given the similarity between its 
CMC properties and those of spiral GMCs? We suggest 
several possibilities. 

1. IC 10 has much more molecular gas than we infer 
from the CO. The discrepancy between IC 10 and large 
galaxies may be entirely explained away by a factor of 
~ 30 increase in the CO-to-Il2 conversion factor. We 
have presented evidence above that BIMA may resolve 
out as much flux as it recovers (leaving us with only a 
factor of ~ 15 discrepancy), but the CMC properties we 
measured above suggest values of Xco that are nearly 
Galactic. Adjusting the data in Figures ^| and El to 
agree with other galaxies requires more than just a large 
reservoir of hidden molecular gas; despite producing rela- 
tively small amounts of CO emission, this molecular gas 
must be associated with star formation. It is unlikely 
that the physical conditions in a hidden reservoir of H2 
could simultaneously be conducive to star formation, in- 
hospitable to CO, and not contribute to the virial masses 
of the GMCs. 

2. The star formation rate of IC 10 is lower than we esti- 
mate from the Ha due to line-blanketing effects in stellar 
atmospheres. These effects would cause us to underpre- 
dict the UV radiation generated per star, leading to an 
overestimate of the star formation rate . Stud ying a sam- 
ple of stars in the SMC. iMassev et alJ (|2005D found that 
O stars in that system (with Z ^ Q.IZq) had effective 
temperatures ^ 4000 K higher than their solar metal- 
licity analogues. For an 05 star {Teff ~ 45,000), this 
results in the SMC star producing ~ 50% more ionizing 
photons. A similar result is obtained using population 
synthesis codes such as STARBURST99 (Figure 78 in 
^cithcrcr et al. 199^ — a shift of an order of magnitude 
in mctallicity for a continuously star forming system re- 
sults in an increase of ~ 50% in the number of ioniz- 
ing photons produced. These adjustments are not large 
enough to make a significant dent in the discrepancy be- 
tween IC 10 and larger galaxies. 

3. If IC 10 has an unexpectedly top heavy IMF 
our inferred star form ation rates may be too high. 
iKennicutt et alJ ()1994|) explores the effect of ad opting 
other (Milky Way) IMFs (such as those of lScalol ^I9m 
or Kroupa et al. ( 1993)) and finds increases of a factor of 
~ 2 — 3 to the star formation rate per imit Ha emission. 
A very top-heavy IMF, by contrast, would have the ef- 
fect of lowering the amount of star formation per unit Ha 
luminosity. The lion's share of ionizing photons are pro- 
duced by stars with M > 10 Mf:, For the Salpeter IMF 
assumed in the Ha calculation fKennic utt et al.l 119941) . 
^ 10% of the mass of stars resides in stars with M > 10 
M0. If IC 10 produced only these stars then the star 
formation efficiency might be a factor of 10 lower than 
the value we calculate here. A very top heavy IMF could 
also explain the unusually high abundance of Wolf Rayet 
stars in IC 10. However, an IMF that is dramatically 
skewed towards high masses would contradict the finding 
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Fig. 13. — Star formation surface density as a func tion of H2 sur- 
face d ensity in IC 10. Global values for large ga laxies |Murgia et alj 
120021) and LMC-size dwarfs ( Lero y et a l."'2 009) are shown for com- 
parison. We convolve the data to a resolution of 1' (275 pc) and 
apply an inclination correction of 0.67. We plot only points in the 
stellar disk of IC 10, where S* > 50 M© pc"^. IC 10 has much 
more star formation per unit molecular gas than most spirals or 
dwarfs. We highlight two other dwarfs similar to IC 10: NGC 1569 
and the SMC. 



bv lHunted ()2001|) that the clusters in IC 10 are consistent 
with a Galactic IMF. Barring such a dramatic IMF, it 
seems unlikely that the offset we observe is only a result 
of a miscalibration of Ha as a star formation tracer. 

4. The star formation rate in IC 10 may have been higher 
in the recent past. If IC 10 has recently undergone a 
period of intense star formation and is now forming stars 
at (relatively) more modest rate, we may be catching 
it at a point its life-cycle during which it is (relatively) 
depleted in molecular gas but still showing the signs of 
a recent star burst. This scenario could also explain the 
very high WR star counts and perhaps the discrepancy 
between the various star formation tracers. In this case, 
IC lO's high star formation efficiency may be temporary, 
an artifact of when we are observing the galaxy. In such a 
case we would expect a large sample of IC 10-like dwarfs 
to average to a position consistent with the other galaxies 
in figures lT^ and ini The higher star formation rate must 
have occurred within the last 10 Myr because high 
mass stars formed during the period of higher SFR must 
still be contributing UV photons that create the Ha flux. 
This is consis tent with the ages of the clusters found by 
IHunted (l20?nh . 4-30 Myr. 

5. Finally, the efficiency of star formation within molec- 
ular clouds may indeed be higher in IC 10 than in the 
Milky Way or other galaxies. This result seems to con- 
tradict the similarities between GMCs in ICIO and those 
in M 31, M 33, and the Milky Way. However, we have 
emphasized the environmental differences between IC 10 
and these systems and these differences may be mani- 
festing themselves in an unexpected way that dramati- 
cally affects the efficiency with which molecular gas forms 
stars. 

6. SUMMARY AND CONCLUSIONS 
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^Hi+Ha (Me pc ^) 

Fig. 14. — The S chmidt Law in IC 10. The points and relation 
of IKennicut^ twMl are shown for comparison. The data have 
been convolved to a resolution of 1' (275 pc) and an inclination 
correction of 0.67 is applied. We plot only points in the stellar 
disk of IC 10, where S« > 50 Mq pc~^. IC 10 shows more star 
formation per unit gas than large galaxies. 



We present a complete survey ^^CO J =(1^0) in 
IC 10 using BIMA. The survey covers all of the opti- 
cal disk of IC 10 and a large part of the extended H I 
structure with a resolution of 14" (70 pc) and a sensitiv- 
ity sufficient to detect clouds with masses greater than 
4 X 10"* Mq. We find structures across the optical disk 
of the galaxy and in the extended structure to the north 
and east of the galaxy. The BIMA finds a total CO lu- 
minosity of 1.0 X 10^ K km s^^ pc^. 

We also present ARO 12m observations of 22 fields in 
IC 10, corresponding to most of the locations in which 
CO emission is detected by the BIMA survey and a num- 
ber of locations of interest with no CO emission. The 
ARO 12m detects CO emission only where BIMA also 
detects CO emission. The ARO 12m finds more emis- 
sion than BIMA or OVRO along the same line of sight. 
This may be evidence for an extended CO component 
surrounding the more compact structures detected by 
the interferometers. Comparing the integrated luminos- 
ity from the ARO 12m to the BIMA survey, we estimate 
that the true CO luminosity of IC 10 is 2.2 x 10^ K km 

— 1 2 

s pc . 

We measure the properties of 14 resolved CO struc- 
tures in IC 1 from high resolution OVRO data by 
IWalter I (l200l . The sizes, line widths, and luminosi- 
ties of these structures resemble thos e of G MCs found 
in similar surveys of M 31, M 33 ijRosolowskv et all 
12001 m.osolowskv fc Lerovl2005D . and Milky Way GMCs 
(jSolomon et al.lll987j) . We conclude that we are observ- 
ing GMCs in IC 10. The virial-mass to luminosity ratio 



in these GMCs is comparable to that observed for spi- 
ral galaxy GMCs and we argue that this implies that 
a Galactic CO-to-H2 conversion factor applies to IC 10. 
We can not constrain the conversion factor in the gas 
resolved out by the interferometers. 

Most of the CO emission detected by the BIMA survey 
comes from lines of sight with T,hi above 10 Mq pc~^ 
{N(H. I) = 1.25 X 10^1 cm-2) and aU of the CO emission 
is very close to such regions. However, only 30% of the 
lines of sight with 'Ehi > 10 Mq pc~^ have associated 
CO emission. This may be because not all high column 
density H I actually corresponds to high volume density 
gas. Indeed, we show that hydrostatic gas pressure pre- 
dicts the CO along a line of sight better than the H I 
column alone. Further, the ratio of molecular to atomic 
gas along a line of sight obeys the same simple relation- 
ship in IC 10 that is seen in large spiral galaxies. 

IC 10 displays significantly more star formation (Ha) 
for a given gas surface density (molecular or total) than 
large spiral galaxies. We suggest several explanations 
for this: a higher CO-to-H2 conversion factor, a different 
IMF or SFR calibration, or a biased (in time) perspec- 
tive. The last explanation is appealing — namely that 
IC lO's unusual star formation efficiency is a timing ef- 
fect. We may be observing a starburst just past its peak, 
so that the signatures of star formation are still present 
but the star-forming gas has already been somewhat de- 
pleted. 
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